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ABSTRACT 

We evaluate the direct detection of extrasolar giant planets with a two-aperture nulling infrared inter- 
ferometer, working at angles 9 < X/2B, and using a new 'ratio-of- two- wavelengths' technique. Simple 
arguments suggest that interferometric detection and characterization should be quite possible for plan- 
ets much closer than the conventional inner working angle, or angular resolution limit. We show that 
the peak signal from a nulling infrared interferometer of baseline (< 40 meters) will often occur 'inside 
the null', and that the signal variations from path-difference fluctuations will cancel to first order in 
the ratio of two wavelengths. Using a new interferometer simulation code, we evaluate the detectability 
of all the known extrasolar planets as observed using this two-color method with the proposed Fourier 
Kelvin Stellar Interferometer (FKSI). In its minimum configuration FK SI uses two 0.5-meter apertures 
on a 12.5-meter baseline, and a ±20° field-of-regard. We predict that ~ 7 known extrasolar planets are 
directly detectable using FKSI, with low-resolution spectroscopy (R ~ 20) being possible in the most 
favorable cases. Spaceborne direct detection of extrasolar giant planets is possible with ~ 12 meter 
baselines, and does not require the much longer baselines provided by formation flying. 

Subject headings: stars: planetary systems — stars: circumstellar matter — techniques: interferometric 
— techniques: high angular resolution — stars: individual (55 Cancri) 



1. INTRODUCTION 

The number of known extrasolar giant planets is now 
over 100 (Marcy et al. 2003), almost all of which 
have been detected indirectly, using the radial velocity 
method. Direct detections of extrasolar giant planets are 
now becoming possible. For example, Charbonneau et al. 
(2002) detected the atmosphere of the 'transiting planet' 
HD 209458b using transit spectroscopy (Seager & Sasselov 
2000). In the case of planets which do not transit, other 
techniques of direct detection are being developed, such as 
visible coronagraphic imaging (Kuchner & Spergel 2003), 
and nulling infrared (IR) interferometry (Bracewell 1978; 
Angel & Woolf 1997). A conventional view is that the 'in- 
ner working angle' (IWA) of a nulling interferometer is the 
angular distance to the first fringe maximum, i.e. X/2B, 
where B is the interferometer baseline and A is wavelength. 
In this Letter, we show that, when the likely IR intensi- 
ties of the extrasolar planets are factored in, nulling IR 
interferometers having B < 40 m, are most sensitive to 
extrasolar giant planets (EGPs) well inside the conven- 
tional IWA, i.e. 'inside the null'. A substantial number 
of the known extrasolar planets are directly detectable us- 
ing interferometry with modest baselines, e.g., the 12.5 m 
baseline of the minimum configuration of FKSI (Danchi et 
al. 2002,2003). 

2. INTERFEROMETRIC DETECTION 'INSIDE THE NULL' ? 

Simple arguments suggest that interferometric detection 
can occur for close-in extrasolar giant planets (CEGPs) not 
fully resolved spatially by the interferometer. First, it is 
well known that a significant source of noise for a space- 
borne nulling interferometer will be photon noise from the 
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leakage of stellar radiation around the 'edges' of the null 
fringe, due to the finite angular radius of the star. If leak- 
age from one stellar radius is a significant noise source, the 
IR signal from the closest CEGPs, at tens of stellar radii, 
will also leak through the null fringe. Second, we note that 
ground-based studies of the CEGPs (Richardson, Deming 
& Seager 2003) are beginning to achieve the sensitivity 
needed to detect the planets in combined light, with no 
spatial resolution whatsoever. 

For a given interferometer baseline, at what angular sep- 
aration does the transmitted IR signal from an extrasolar 
planet peak? For decreasing angular separation from the 
star, the fringe transmission decreases, but the planet's 
thermal equilibrium temperature in the stellar radiation 
field increases, so the signal transmitted through the null 
fringe can remain significant. Distant planets will be ob- 
served with fringe transmission of unity, but they will be 
colder. Figure 1 shows the product of the Planck function 
at a wavelength of 5 /*m, times the fringe transmission 
(sin 2 (7r6»B/A)) for baselines of 8, 12, 20 and 40 m. The ex- 
trasolar system was assumed to be at 10 pc distance, and 
the planet was assumed to emit as a blackbody in thermal 
equilibrium with a Bond albedo of 0.4. These 'interfer- 
ometer contribution functions' were normalized to unity 
for the peak signal at the 40-m baseline. The asterisks on 
each curve mark the nominal IWA at X/2B. 

Figure 1 shows that, even for the longest (40-m) base- 
line, the peak signal occurs inside the nominal IWA at 
X/2B. For the shorter baselines the peak occurs quite far 
inside the IWA. For example, the 12-m peak signal occurs 
for planets near 0.1 AU, a factor of 4 below X/2B. For all 
baselines, the decrease in signal at the greatest distances 
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occurs because of the lower planetary temperatures. Al- 
though the peak signal for the 12-m baseline occurs near 
0.1 AU, this signal is reduced by about an order of magni- 
tude from the 40-m case. Figure 1 suggests that the sen- 
sitivity of a 12-m interferometer working 'inside the null' 
may be sufficient to allow detection and characterization 
of some of the known extrasolar giant planets. 

3. A TWO-COLOR METHOD 

Using a simple analytical model of the nulling interfer- 
ometer we demonstrate that the ratio of the intensity at 
the output of the nuller at two wavelengths is insensitive 
to the residual pathlength fluctuations, provided they are 
small compared to the wavelength. A nulling interferom- 
eter operates like a conventional stellar interferometer ex- 
cept that an achromatic tt phase shift is applied to the 
beam coming from one element of the interferometer and 
that the beams are symmetrically and achromatically com- 
bined. The equation for the normalized intensity, N(Xi), 
is: 



JV(Ai) = l/2 [l-|V(Ai)| cos «£(Ai)] 



(1) 



where Ai is the first wavelength, |V(Ai)| is the modu- 
lus of the visibility of the source, and </>(Ai) = 2-Ka^j\\ 
is the rms phase error caused by residual pathlength fluc- 
tuations, <jn in the system. For a point source and no 
pathlength fluctuations, Eq. 1 reduces to the usual 9 2 null 
of the classical Bracewell interferometer. A similar equa- 
tion holds for the output of the same nuller at a second 
wavelength, A2. Using these assumptions, 



(2) 



where 8 sta r is the angular diameter of the star and B is 
the baseline length, and similarly expanding the term for 
the phase fluctuations: 

cos (/.(Ai) w 1 - <?>(Ai)72 = 1 - 2tt 2 (oat/Ai) 2 (3) 

Substituting these equations into (1) and using similar 
equations for the second wavelength, A2 , it is easy to show 
that: 



iV(Ai) 
A(A 2 ) 



(4) 



Hence the stellar leakage and pathlength fluctuations 
cancel out. 

This result is valid for the case when the interferometer 
is not rotating about the line of sight or when the residual 
pathlength fluctuations occur at much higher frequencies 
than the rotation frequency about the line of sight. The 
cancellation occurs in the leakage signal variations which 
accompany instability in the null fringe (not in the photon 
noise), and which would otherwise overwhelm the plane- 
tary signal. This technique is expected to work best when 
the two wavelengths Ai and A2 are reasonably close to- 
gether, avoiding higher order effects. 

4. AN INTERFEROMETER SIMULATION CODE 

We have written a simulation code to compute the sig- 
nal from the known extrasolar giant planets, as observed 



by a rotating 2-aperture nulling interferometer, following 
the principle of Bracewell (1978), and treating the out- 
put signal as a ratio, exploiting the rationale given above. 
The code uses the interferometer parameters for FKSI as 
listed in Table 1. The modeled fringe pattern of the in- 
terferometer includes path difference errors, based on the 
current FKSI error budget. Shot-noisc is computed from 
the total signal, including the stellar leakage, the extra- 
solar zodiacal background, and the instrument thermal 
background. Dark current and read noise from the de- 
tection system are also included. The stellar spectrum 
is assumed to be a blackbody, but wavelength-dependent 
limb darkening is imposed on the disk, based on the solar 
observations of Pierce et al. (1950). The stellar tem- 
perature and radius are estimated from the spectral type. 
The planet is modeled as a blackbody spectrum with su- 
perposed molecular band spectral structure. The planet's 
blackbody equilibrium temperature is computed assuming 
a Bond albedo of 0.4; the planetary radius is taken to be 
35% greater than Jupiter, based on HD 209458b (Brown et 
al. 2001). The planetary spectral structure was included 
by interpolating in the 'cloudless' sequence computed by 
Sudarsky, Burrows & Hubeny (2003). The planetary spec- 
trum is currently assumed to be independent of orbital 
phase and inclination to the line of sight. Orbital motion 
of the planet is included, using orbital elements from the 
Doppler observations, and a nominal inclination of 45°. 
The zodiacal background is calculated including contri- 
butions from both scattered light and thermal radiation 
(Kuchncr 2002). Since younger stars will have more mas- 
sive zodiacal disks, we scaled the mass of the disk as the 
-1.76 power of the stellar age (Spangler et al. 2001). Ages 
for stars hosting extrasolar planets were taken from Laws 
et al. (2003). The inclination of the zodiacal disk is also 
included (nominally 45° to the plane of the sky), since the 
asymmetry from the inclined disk produces a significant 
signal as the interferometer rotates. 

5. EXAMPLE OF AN FKSI DETECTION 

We have simulated observations by FKSI of all known 
extrasolar planetary systems using the code. An exam- 
ple is shown in Fig. 2, which illustrates the detection of 
55 Cancri b. Figure 2a shows the number of photons de- 
tected versus wavelength, at the FKSI spectral resolving 
power (X/8X = 20), in one 300-s integration. The calcu- 
lated wavelength range in Fig. 2 extends longward of the 
nominal wavelength limits of FKSI. The dominant source 
of photons is stellar leakage, with planetary radiation sec- 
ond. The zodiacal radiation falls below the planetary in- 
tensity for this old, and dust-poor system (Jayawardhana 
et al. 2002). The dominant source of noise depends on 
wavelength, but is due to stellar leakage shot noise at the 
shortest wavelengths, thermal background radiation noise 
at the longest wavelength (8 fim), and dark current and 
read noise at intermediate wavelengths. Since the antici- 
pated total integration time during an observational 'cam- 
paign' for a given planetary system can be many days, 
planets whose photon counts fall well below the noise level 
in a single 300-s integration will nevertheless be detectable. 
But for 55 Cancri b, the planetary signal photon counts 
are comparable to the noise photons in a single 300-s inte- 
gration, for wavelengths > 3 fim. This planet is strongly 
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detectable, and more extensive spectral information could 
also be extracted. 

In the original Bracewell concept, the rotation of the in- 
terferometer modulates the planetary signal, and we have 
included this process in the simulation code. Figure 2b 
(lower panel) shows a power spectrum from a signal time 
series of 55 Cancri b. As per our two-color method, the 
Fourier-transformed signal in this case was the intensity 
integrated over wavelengths 6 ^ A ^ 3 /um, ratioed to 
the intensity integrated over A < 3 /xm. The denomina- 
tor in this ratio is dominated by the stellar leakage at the 
shortest wavelengths, and contains minimal planetary sig- 
nal. The planet peaks at the longer wavelengths in the 
numerator. The greater null depth at long wavelengths 
enhances this separation of planet and stellar signals. 

The interferometer was rotated slowly, about 15 hours 
per rotation, and the simulated campaign lasted about 35 
days (50 rotations). The zodiacal signal appears in the 
Fig. 2b power spectrum at twice the interferometer rota- 
tion frequency, i.e. at 37 microHz (/iHz), and the plane- 
tary signal also contributes at this frequency. Because of 
the planet's orbital motion, the planetary signal appears 
primarily at a different frequency (see caption), slightly 
displaced, and resolved from the zodiacal signal. More- 
over, because the planet is inside the null, no overtones 
are caused by transmission through higher-order fringes. 
Note also that, as the signal from the planet goes in- and 
out-of-phase with the zodiacal signal, the envelope of this 
modulation is seen in the low- frequency region of the power 
spectrum as a peak at 1.6 /xHz, twice the orbital frequency 
of the planet. 

6. DISCUSSION 

Our simulations show that a significant number of 
known extrasolar planets will be detectable using FKSI, 
and that spectral information can be derived in several 
of the most favorable cases. These simulated detections 
are robust in the sense that the few most favorable plan- 
etary systems are detectable with almost any reasonable 
instrumental configuration (baseline, aperture size, etc.) 
or choice of long-wavelength cutoff. However, detections 
for the fainter systems depend significantly on the instru- 
mental parameters and wavelength range, as well as on the 
properties of the detectors (dark current, read noise). 

Figure 3 illustrates the distances and orbital semi-major 
axes for many of the known planetary systems, with those 
systems detectable by FKSI (using Table 1 parameters) 
plotted as filled symbols. The ±20° ficld-of-regard {FOR) 
holds seven detectable planets, and an additional six plan- 
ets are within a ±40° FOR. The conventional IWA resolu- 



tions for baselines of 40, 20, 12, and 8 m are overplotted. 
It can be seen that the detectable planets are essentially 
the CEGP systems closest to earth, with essentially no 
dependence on the angular resolution limits (IWA lines). 

Our results show that planets can be detected much 
closer to stars than expected based on the concept of the 
inner working angle, computed from the nominal resolu- 
tion of the interferometer, i.e., X/2B. This resolution es- 
timate is essentially the same as the Rayleigh criterion for 
conventional telescopes, which is the angular separation of 
two stars of equal intensity in which one star is placed at 
the first zero of the Airy pattern of the second star. The 
Rayleigh criterion is well known to be a very conservative 
estimate of resolution, and sources can be resolved that are 
substantially closer than this, using super-resolution tech- 
niques such as the pixon method (Pina & Puetter 1993). 
Similar considerations hold for interferometers as CLEAN, 
MEM, and other methods have provided images with ef- 
fective angular resolutions much better than nominally ex- 
pected (Cobb & Fix 1987). 

For a wavelength of 5 fim, the resolution of the FKSI 
is X/2B « 41 milli-arcsec (mas). But our simulations 
demonstrate an 'effective resolution' of approximately 1 
mas. Our results are consistent with the working resolu- 
tion of the interferometer being determined by the ratio 
of the rms pathlength fluctuation to the baseline, <jn/B, 
which for the parameters in this paper, is 0.2 mas. This 
is not as surprising as it may seem at first, since the path- 
length stability requirement to achieve a 10~ 4 null or bet- 
ter is essentially that of an interferometer with extremely 
high phase stability. 

The work presented here has important implications for 
the Terrestrial Planet Finder (TPF) mission (Bcichman, 
Woolf, & Lindensmith 1999), because angular resolution is 
often emphasized over sensitivity. The desire for a nulling 
interferometer using free flyer telescopes is based on the 
resolution needed to search more than 150 F, G, and K 
stars for earth-like planets in the habitable zone (Luninc 
et al. 2003). Given that the actual IWA can be signifi- 
cantly smaller than previously thought, means that it may 
be possible to achieve the basic goals of TPF with a struc- 
turally connected interferometer having a modest baseline 
in the range of 20 to 30 m. 

7. ACKNOWLEDGEMENTS 

We thank Drs. R. Allen, D. Bcnford, D. Gczari, D. Lei- 
sawitz, J. Monnier, M. Mumma, L. Mundy, C. Noecker, 
and W. Traub, for their contributions to the FKSI mis- 
sion concept. 



REFERENCES 



Angel, J. R. P. & Woolf, N. J. 1997, ApJ, 475, 373 

Beichman, C, Woolf, N. J., & Lindensmith, C. 1999 (eds.), JPL 

Publication 99-3. 
Bracewell, R. N. 1978, Nature, 274, 780 

Brown, T. M., Charbonneau, D., Gilliland, R. L., Noyes, R. W., & 

Burrows, A. 2001, ApJ, 552, 699 
Charbonneau, D., Brown, T. M., Noyes, R. W., & Gilliland, R. L. 

2002, ApJ, 568, 377 
Cobb, M. L., & Fix, J. D. 1987, ApJ, 315, 325 
Danchi, W. C, ct al. 2002, AAS Meeting 201, # 151.09 
Danchi, W. C. et al. 2003, in Towards Other Earths , ESA 

Publication SP-539, in press. 



Jayawardhana, R., Holland, W. S., Kalas, P., Greaves, J. S., Dent, 

W. R. F., Wyatt, M. C, & Marcy, G. W. 2002, ApJ, 570, L93 
Kuchner, M. J. 2002, IDL software program 'zodipic', http:// cfa- 

www.harvard.edu/~mkuchncr 
Kuchner, M. J., and Spergel, D. 2003, in Scientific Frontiers in 

Research on Extrasolar Planets, ASP Conference Series vol. 294 

(D. Doming & S. Seager, eds.), p. 603 
Laws, C, Gonzalez, G., Walker, K. M., Tyagi, S., Dodsworth, J., 

Snider, K., & Suntzeff, N. B. 2003, AJ, 125, 2664 
Lunine, J. I. et al. 2003, TPF Science Working Group, unpublished. 
Marcy, G. W., Butler, R. P., Fischer, D. A., & Vogt, S. S. 2003, 

in Scientific Frontiers in Research on Extrasolar Planets, ASP 

Conference Series vol. 294 (D. Deming k, S. Seager, eds.), p. 1 



4 



Danchi et al. 



Pierce, A.K., McMath, R. R., Goldberg, L., & Mohlcr, O. C. 1950, 

ApJ, 112, 289 
Pina, R. K., & Puetter, R. C. 1993, PASP, 105, 630 
Richardson, L. J., Doming, D., & Scager, S. 2003, ApJ, in press, 

astro-ph/0307297 



Seager, S. & Sasselov, D. D. 2000, ApJ, 537, 916 

Spangler, C, Sargent, A. I., Silverstone, M. D., Becklin, E. E., 

Zuckerman, B. 2001, ApJ, 555, 932 
Sudarsky, D., Burrows, A., & Hubeny, I. 2003, ApJ, 588, 1121 



Extrasolar Giant Planets Using FKSI 



5 



Table 1 

Nominal Parameters for FKSI 



Parameter 



Value 



Baseline 
Aperture 
Field of regard 
Efficiency 
Spectral resolution 
Pathlcngth stability 
Wavelength range 
Optics temperature 
Detector temperature 
Dark current 
Read noise 



12.5 meters 

2 x 0.5 meters diameter 
±20° 

0.05 (electrons out / photons in) 
\/5\ = 20 

15 nm rms; 10-sec time constant 

1 ^m ^ A ^ 6 fim 

63K 

35K 

0.2 e^/sec 
8 e" 
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planet distance (AU) 



Fig. 1. — Log of the peak planetary signal seen at A = 5 fim by a nulling interferometer of (upper to lower curves) baseline 40, 20, 12, and 8 
meters. The asterisks mark the nominal inner working angles at X/2B. The curves were normalized to unity for the peak signal at the 40-m 
baseline. The IR intensity of the planet was computed from blackbody equilibrium, with Bond albedo =0.4, and distance 10 pc. 
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Fig. 2. — Simulated observations of the innermost planet in the 55 Cancri system (a = 0.11 AU, P = 14.65 days) using FKSI. The upper 
panel (Figure 2a) shows the number of photons detected at a single interferometer rotation angle (where the planet signal is maximum), in 
a 300 second integration. The stellar leakage is the upper solid line, and the detected planetary spectrum is the lower solid line. The dashed 
line is the zodiacal component, and the dotted line is the total noise. The lower panel (Figure 2b) shows the power spectrum from a simulated 
35-day observing campaign (50 rotations of the interferometer). The marks on the frequency axis indicate 2i>i (solid mark), and 2\vi — u p \ 
(dashed mark), where Vi is the rotation frequency of the interferometer, and v v = 1/P, where P is the orbital period of the planet. 
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Fig. 3. — Distance (pc) and orbital radius (semi-major axes, in AU) of known extrasolar planets (having d ^ 50 pc, and r ^ 4 AU, and 
within a ±40° field of regard). The planets nominally detectable by FKSI are plotted with filled symbols; circles indicate planets within a 
±20° field of regard, and diamonds indicate additional detectable planets if the field of regard is extended to ±40° . The lines correspond to 
the first fringe maximum for nulling interferometers of (top to bottom) baselines 8, 12, 20, and 40 meters, at A = 5 jim. 



